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ABSTRACT 

Understanding infrared (IR) luminosity is fundamental to understanding the cosmic star for¬ 
mation history and AGN evolution, since their most intense stages are often obscured by dust. 
Japanese infrared satellite, AKARI, provided unique data sets to probe this both at low and high 
redshifts. The AKARI performed all sky survey in 6 IR bands (9, 18, 65, 90, 140, and 160/mr) 
with 3-10 times better sensitivity than IRAS, covering the crucial far-IR wavelengths across the 
peak of the dust emission. Combined with a better spatial resolution, AKARI can much more 
precisely measure the total infrared luminosity ( Ltir ) of individual galaxies, and thus, the total 
infrared luminosity density of the local Universe. In the AKARI NEP deep field, we construct 
restframe 8/im, 12/rm, and total infrared (TIR) luminosity functions (LFs) at 0.15< z <2.2 using 
4128 infrared sources. A continuous filter coverage in the mid-IR wavelength (2.4, 3.2, 4.1, 7, 9, 11, 
15, 18, and 24/um) by the AKARI satellite allows us to estimate restframe 8/im and 12/rm lumi¬ 
nosities without using a large extrapolation based on a SED fit, which was the largest uncertainty 
in previous work. By combining these two results, we reveal dust-hidden cosmic star formation 
history and AGN evolution from z= 0 to z=2. 2, all probed by the AKARI satellite. 
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1. INTRODUCTION 

Revealing the cosmic star formation history is one of 
the major goals of the observational astronomy. How¬ 
ever, UV/optical estimation only provides us with a 
lower limit of the star formation rate (SFR) due to 
the obscuration by dust. A straightforward way to 
overcome this problem is to observe in the infrared, 
which can capture the star formation activity invis¬ 
ible in the UV. The superb sensitivities of recently 
launched Spitzer and AKARI satellites can revolution¬ 
ize the held. 

In the local Universe, often used IR LFs are from the 
IRAS (e.g., Sanders et al. 2003; Goto et al. 2011a) from 
1980s, with only several hundred galaxies. In addition, 
bolometric infrared luminosities (L/r, 8-ioooAtm) of lo¬ 
cal galaxies were estimated using equation in Perault 
(1987), which was a simple polynomial, obtained as¬ 


suming a simple blackbody and dust emissivity. Fur¬ 
thermore, the reddest filter of IRAS was 100/im, which 
did not span the peak of the dust emission for most 
galaxies, leaving a great deal of uncertainty. Using 
deeper AKARI all sky survey data that cover up to 
160/im, we aim to measure local L/r, and thereby the 
IR LF more accurately. 

At higher redshifts, most of the Spitzer work relied 
on a large extrapolation from 24/mr flux to estimate the 
8, 12/im or total infrared (TIR) luminosity, due to the 
limited number of mid-IR Liters. AKARI has continu¬ 
ous Liter coverage across the mid-IR wavelengths, thus, 
allows us to estimate mid-IR luminosity without using 
a large fc-correction based on the SED models, elim¬ 
inating the largest uncertainty in previous work. By 
taking advantage of this, we present the restframe 8, 
12/rm TIR LFs, and thereby the cosmic star formation 


http://pkas.kas.org 


1 


2 


Goto 



Figure 1: Emission line ratios used to select AGNs 
from the AKARI all sky sample. The contour shows 
distribution of all galaxies in the SDSS with r < 17.77 
(regardless of IR detection). The dotted line is the 
criterion between starbursts and AGNs described in 
Kewley et al. (2001). The dashed line is the criterion 
by Kauffmann et al. (2003). Galaxies with line ra¬ 
tios higher than the dotted line are regarded as AGNs. 
Galaxies below the dashed line are regarded as star¬ 
forming. Galaxies between the dashed and dotted lines 
are regarded as composites. The blue and red dots are 
for ULIRGs, LIRGs, respectively. The green squares 
are Seyfert 1 galaxies identified by visual inspection of 
optical spectra. 


history derived from these using the AKARI NEP-Deep 
data. 

2. AKARI All Sky Survey: low-;; Universe 

AKARI performed an all-sky survey in two mid- 
infrared bands (centered on 9 and 18 /mi) and in four 
far-infrared bands (65, 90, 140, and 160^tm). In addi¬ 
tion to the much improved sensitivity and spatial reso¬ 
lution over its precursor (the IRAS all-sky survey), the 
presence of 140 and 160/rm bands is crucial to measure 
the peak of the dust emission in the FIR wavelength, 
and thus the Ljr of galaxies. We have cross-correlated 
the AKARI FIS bright source catalog with the SDSS 
DR7 galaxy catalog, obtaining 2357 spectroscopic red- 



Figure 2: Fractions of AGN and composite galaxies as 
a function of L/_r. AGN are classified using Kewley et 
al. (2001) among galaxies with all 4 lines measured. 
Composite galaxies include those classified as AGN us¬ 
ing Kauffmann et al. (2003). 
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Figure 3: Infrared luminosity function of AKARI- 
SDSS galaxies. The Lj^ is measured using the AKARI 
9, 18, 65, 90, 140 and 160/xm fluxes through an SED 
fit. Errors are computed using 150 Monte Carlo simula¬ 
tions, added to a Poisson error. The dotted lines show 
the best-fit double-power law. The green dotted lines 
show IR LF at 2=0.0082 by Goto et al. (2011a). The 
dashed-dotted lines are higher redshift results from the 
AKARI NEP deep field (Goto et al. 2010a,b). 
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Figure 4: (left) Restframe 8^m LFs. The blue diamonds, purple triangles, red squares, and orange crosses show 
the 8/rm LFs at 0.38 < 2 < 0.58,0.65 < z < 0.90,1.1 < 2 < 1.4, and 1.8 < 2 < 2.2, respectively. The dotted lines 
show analytical fits with a double-power law. Vertical arrows show the 8/rm luminosity corresponding to the flux 
limit at the central redshift in each redshift bin. Overplotted are Babbedge et al. (2006) in the pink dash-dotted 
lines, Caputi et al. (2007) in the cyan dash-dotted lines, and Huang et al. (2007) in the green dash-dotted lines. 
AGNs are excluded from the sample, (middle) Restframe 12^m LFs. The blue diamonds, purple triangles, and red 
squares show the 12/rm LFs at 0.15 < z < 0.35,0.38 < z < 0.62, and 0.84 < z < 1.16, respectively. Overplotted are 
Perez-Gonzalez et al. (2005) at z=0.3,0.5 and 0.9 in the cyan dash-dotted lines, and Rush et al. (1993) at z=0 in 
the green dash-dotted lines, (right) TIR LFs. The redshift bins used are 0.2< z <0.5, 0.5< z <0.8, 0.8< z <1.2, 
and 1.2< z <1.6. Overplotted are Le Floc’h et al. (2005), Magnelli et al. (2009), Huynh et al. (2007), and Sanders 
et al. (2003). 
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Figure 5: Evolution of TIR luminosity density based on TIR LFs (red circles), 8/xm LFs (stars), and 12/im LFs 
(filled triangles). The blue open squares and orange filled squares are for LIRG and ULIRGs only, also based on 
our Ltir LFs. Overplotted dot-dashed lines are estimates from the literature: Le Floc’h et al. (2005), Magnelli 
et al. (2009), Perez-Gonzalez et al. (2005), Caputi et al. (2007), and Babbedge et al. (2006) are in cyan, yellow, 
green, navy, and pink, respectively. The purple dash-dotted line shows UV estimate by Schiminovich et al. (2005). 
The pink dashed line shows the total estimate of IR (TIR LF) and UV. 
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Figure 6: Evolution of TIR luminosity density by 
AGN. Results from the AKARI all sky survey is plot¬ 
ted with stars at £=0.0082. The red, blue and orange 
points show IR luminosity density from all AGN, from 
LIRG AGN only, and from ULIRG AGN only. Higher 
redshift results are from the AKARI NEP deep field 
(Goto et al. 2010a), with contribution from star form¬ 
ing galaxies removed. Brown triangles are com¬ 

puted from the 8^m LFs (Goto et al. 2010a). 

shifts. 

It is fundamental to separate IR contribution from 
two different physics; the star-formation and AGN ac¬ 
tivity. In Fig. [TJ we use [NII]/Ha vs [OIII]/H/3 
line ratios to classify galaxies into AGN or SFG (star¬ 
forming galaxies). It is interesting that majority of 
(U)LIRGs are aligned along the AGN branch of the 
diagram, implying the AGN fraction is high among 
(U)LIRGs. This is more clearly seen in Fig. [2J where 
we plot fractions of AGN as a function of Lip. This re¬ 
sults agree with previous AGN fraction estimates (Goto 
et al. 2005; Yuan et al. 2010). Improvement in this 
work is that due to much larger statistics, we were able 
to show fractions of AGN in much finer luminosity bins, 
more accurately quantifying the increase. Especially, a 
sudden increase of f agn at log Lip >11.3 is notable 
due to the increased statistics in this work. 

For these galaxies, we estimated total IR luminosi¬ 
ties (Lip) by fitting the AKARI photometry with SED 
templates. We used the LePhare cod^j] to fit the in¬ 
frared part (>7/mr) of the SED. We fit our AKARI 

Http: //www. cfht.hawaii.edu/~arnouts/lephare.html 


Table 1: Summary of evolution of infrared luminosity. 



oc(l+z ) 4 - 1±0 ' 4 

(ULIRG) 

oc(l+z ) 100±0 - 5 

nfFG(LiRG) 

oc(l+z ) 6 - 5±0 - 5 


oc(l+z ) 4 - 1±0 - 5 

nf£ N (ULIRG) 

oc(l+z ) 8 - 7±0 - 6 

nf° N (LIRG) 

oc(l+z ) 5 - 4±0 - 5 


FIR photometry with the SED templates from Chary 
& Elbaz (2001; CHEL hereafter), which showed most 
promising results among SED models tested by Goto 
et al. (2011a). 

With accurately measured Lip, we are ready to con¬ 
struct IR LFs. Since our sample is flux-limited at 
r = 17.7 and Sgo^m = 0.7 Jy, we need to correct for 
a volume effect to compute LFs. We used the l/V max 
method. We estimated errors on the LFs with 150 
Monte Carlo simulations, added to a Poisson error. 

In Fig. [3] we show infrared LF of the AKARI-SDSS 
galaxies. The median redshift of our sample galaxies is 
£=0.031. 

Once we measured the LF, we can estimate the 
total infrared luminosity density by integrating the 
LF, weighted by the luminosity. We used the best- 
fit double-power law to integrate outside the luminos¬ 
ity range in which we have data, to obtain estimates 
of the total infrared luminosity density, LIir. Note 
that outside of the luminosity range we have data 
(Lip > 10 12 5 Lq or Lip < 10 9 - 8 L Q ), the LFs are 
merely an extrapolation and thus uncertain. 

The resulting total luminosity density is Llip= 
(3.8j^; 8 ) x 10 s LqMpc - 3 . Errors are estimated by vary¬ 
ing the fit within la of uncertainty in LFs. Out of Llip, 
1.1A0.1% is produced by LIRG (Lip > 10 11 Lq), and 
only 0.03±0.01% is by ULIRG (Lip > 1O 12 L 0 ). Al¬ 
though these fractions are larger than £=0.0081 (Goto 
et al. 2011a), still a small fraction of flip is produced 
by luminous infrared galaxies at £=0.031, in contrast 
to high-redshift Universe. 

3. AKARI NEP Deep Field: high-£ Universe 

The AKARI has observed the NEP deep field (0.4 deg 2 ) 
in 9 filters (N2, 7V3, IV4,57, S9W, 511, L15, L18W and 
L24) to the depths of 14.2, 11.0, 8.0, 48, 58, 71, 117, 
121 and 275/xJy (5cr; Wada et al. 2008). This re¬ 
gion is also observed in BVRi'z' (Subaru), u' (CFHT), 
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FUV, NUV (GALEX), and J,Ks (KPN02m), with 
which we computed photo-z with ^^=0.043. Objects 
which are better fit with a QSO template are removed 
from the analysis. We used a total of 4128 IR sources 
down to 100 /iJy in the L18 filter. We compute LFs 
using the 1/Vmax method. Data are used to 5 a with 
completeness correction. Errors of the LFs are from 
1000 realization of Monte Carlo simulation. 

3.1. 812 m LF 

Monochromatic 8/im luminosity (Lg^m) is known to 
correlate well with the TIR luminosity (Babbedge et 
al. 2006; Huang et al. 2007; Goto et al. 2011b), espe¬ 
cially for star-forming galaxies because the rest-frame 
8/im flux is dominated by prominent PAH features such 
as at 6.2, 7.7 and 8.6 /cm. The left panel of Fig. |4]shows 
a strong evolution of 8/im LFs. Overplotted previous 
work had to rely on SED models to estimate Lg/im from 
the Spitzer i n the MIR wavelengths where SED 

modeling is difficult due to the complicated PAH emis¬ 
sions. Here, AKARI’s mid-IR bands are advantageous 
in directly observing redshifted restframe 8/rm flux in 
one of the AKARI’s filters, leading to more reliable 
measurement of 8-/rm LFs without uncertainty from 
the SED modeling. 

3.2. 12/im LF 

12/rm luminosity (Li 2 ^m) represents mid-IR contin¬ 
uum, and known to correlate closely with TIR lumi¬ 
nosity (Perez-Gonzalez et al. 2005). The middle panel 
of Fig. [4] shows a strong evoltuion of 12/rm LFs. Here 
the agreement with previous work is better because (i) 
12/rm continuum is easier to be modeled, and (ii) the 
Spitzer also captures restframe 12/xm in S , 24 M m at 2 =1. 

3.3. TIR LF 

Lastly, we show the TIR LFs in the right panel of Fig. 
|4j We used Lagache et al. (2003)’s SED templates to 
fit the photometry using the AKARI bands at >6/mi 
{S7,S9W,S11,L15,L18W and L24). The TIR LFs 
show a strong evolution compared to local LFs. At 
0.25 < 2 < 1.3, L RIR evolves as oc (1 + z) 4 ' 1±a4 . 

4. Cosmic star formation history 

We fit LFs in Fig. [4] with a double-power law, then 
integrate to estimate total infrared luminosity density 


at various z. The restframe 8 and 12/iin LFs are con¬ 
verted to Ltir using Perez-Gonzalez et al. (2005) and 
Caputi et al. (2007) before integration. The resulting 
evolution of the TIR density is shown in Fig. [5] The 
right axis shows the star formation density assuming 
Kennicutt (1998). We obtain ( 2 ) oc (l + z) 41±0 - 4 . 
Comparison to fluv using Schiminovich et al. (2005) 
suggests that Qtir explains 70% of fltotai at 2=0.25, 
and that by 2=1.3, 90% of the cosmic SFD is explained 
by the infrared. This implies that flriR provides good 
approximation of the Qtot.ai at 2 > 1. 

In Fig. [5] we also show the contributions to LLtir 
from LIRGs and ULIRGs. From 2=0.35 to 2=1.4, Hjr 
by LIRGs increases by a factor of ~1.6, and Hir by 
ULIRGs increases by a factor of ~10. More details are 
in Goto et al. (2010a). 

5. Cosmic AGN accreation history 

We have separated the from flf R N . We can also 

investigate flf R N . 

In Fig. [Gj we show the evolution of flf R N , which 
shows a strong evolution with increasing redshift. At a 
first glance, both f lf R N and Hj R g show rapid evolu¬ 
tion, suggesting that the correlation between star for¬ 
mation and black hole accretion rate continues to hold 
at higher redsliifts, i.e., galaxies and black holes seem 
to be evolving hand in hand. When we fit the evo¬ 
lution with (l+z) T , we found £lf R N oc(l+z) 4 - 1±0 - 5 . 
A caveat, however, is that &f R N estimated in this 
work is likely to include IR emission from host galax¬ 
ies of AGN, although in optical the AGN compo¬ 
nent dominates. Therefore, the final conclusion must 
be drawn from a multi-component fit based on bet¬ 
ter sampling in FIR by Herschel or SPICA, to sepa¬ 
rate AGN/SFG contribution to Ljr. The contribu¬ 
tion by ULIRGs quickly increases toward higher red- 
shift; By 2=1.5, it exceeds that from LIRGs. In¬ 
deed, we found flf R N (ULIRG) oc(l+z) 8 - 7±0 - 6 and 
flf GN (LIRG) oc(l+z) 5 - 4±0 - 5 . Summary of the evolu¬ 
tion of Hir is presentedtj in Table [1] (See Goto et al. 
2011b for more details). 
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